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1.  INTRODUCTION 

A  key  problem  in  the  design  of  communication  networks  is  the  efficient 
sharing  of  a  common  transmission  channel,  (such  as  a  satellite  link,  a  ground 
radio  channel,  a  computer  bus,  a  coaxial  cable,  or  an  optical  fibre)  among  a 
large  population  of  network  users.  This  problem  is  referred  to  as  the  multiple 
access  problem,  since  many  independent  users  share,  and,  thus,  access  a  common 
channel  for  transmission  of  information.  The  solution  to  the  multiple-access 
problem  must  incorporate  a  distributed  control  scheme,  termed  multiple-access 
algorithm,  for  allocating  the  channel  resources  among  the  network  users. 

The  design  and  performance  of  multiple-access  algorithms  are  highly 
dependent  on  the  nature  of  the  users.  When  a  channel  is  to  support  large 
numbers  of  bursty  (low  duty-cycle)  users,  random  multiple-access  algorithms 
( RMAAs )  become  more  efficient  than  deterministic  algorithms.  This  has  been 
early  recognized  by  the  researchers  in  the  field,  and  a  plethora  of  RMAAs 
have  been  proposed  during  the  past  fifteen  years  [1,2] . 

The  key  performance  measures  of  a  RMAA  are  its  throughput  and  delay 
characteristics.  The  evaluation  of  such  characteristics  has  been  the  subject 
of  numerous  studies.  In  most  cases,  a  Markovian  model  is  employed,  and  the 
existence  of  steady  state  of  the  random-access  system  is  related  to  the 
ergodicity  of  an  underlying  Markov  process.  Depending  on  the  complexity  of 
the  state  space  of  such  a  process,  this  formulation  usually  gives  sufficient 
information  on  the  maximum  input  traffic  rate  that  an  algorithm  can  maintain. 
However,  the  evaluation  of  the  delay  charactr r i st ics  is  a  much  harder  problem, 
since  they  are  intimately  ltnerwoven  with  the  dynamical  behavior  of  the 
algorithm's  scheduling  mechanisms.  Due  to  this  fact,  it  is  not  surprising  that 
results  concerning  the  delay  character  ist.  i*'.,  are  limited,  anil  are  obtained  afte 
a  rather  intricate  and  iiffieult  m  ilysis,  which  is  usually  matched  to  the 
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peculiarities  of  the  specific  algorithm  at  hand. 

In  this  paper,  we  show  how  the  delay  analysis  of  RMAAs  can  be  unified  and 
simplified,  by  the  use  of  some  known  results  from  the  theory  of  regenerative 
processes,  and  the  theory  of  infinite  dimensional  systems  of  linear  equations. 
After  outlining  the  method  in  section  2,  we  demonstrate  its  wide  applicability 
and  relative  simplicity,  by  applying  it,  in  sections  3  and  4,  to  three  algorithms 
that  represent  different  classes  of  RMAAs,  namely: 

1)  the  Controlled  ALOHA  algorithm  ( "ALOHA- type"  class)  [6J 

2)  the  "0.487"  algorithm  ("full  sensing-blocked  access"  class)  [7,8] 

3)  the  n-ary  stack  algorithm  ("limited  sensing-free  access"  class) 

[15,16,17] 

flor  the  above  algorithms,  we  obtain  explicit  results  on  the  induced  mean 
delay,  for  the  Poisson  infinite-user  population  model.  The  higher  moments  of 
the  delay,  for  the  Poisson  as  well  as  for  an  arbitrary  memoryless  input  stream, 
can  be  computed  using  the  same  method. 

2 .  THE  METHOD 

In  random-access  systems,  as  in  virtually  every  queueing  system,  many  of 
the  involved  stochastic  processes  are  regenerative.  A  regenerative  process  is 
a  process  consisting  of  a  sequence  of  regeneration  cycles  such  that  the 
probabilistic  structure  of  the  sample  function  of  the  process  during  such  a 
cycle  is  the  same  for  every  cycle,  and  independent  of  that  of  previous  or 
future  cycles  [5].  The  beginnings  of  such  cycles  are  referred  to  as  regeneration 
times ,  and  form  a  renewal  process.  A  discrete-time  process  xXn)n>;i  is 
said  to  be  regenerative  with  respect  to  the  renewal  process  if 

the  process  lx^  +n)n  ^  is  a  probabilistic  replica  of  the  process  [Xn;ni1, 


for  every  i,  i=l,2,... 


The  regenerative  pr>  >babi  1 1  st  ic  structure  of  a  regenerative  process  makes 


it  possible  to  express  its  asymptotic  behavior  in  terms  of  quantities  that 
refer  only  to  one  cycle  of  the  process*  Tliis  is  made  precise  by  the  following 
elegant  and  powerful  result,  which  will  be  referred  to  as  the  regeneration 
theorem.  [3,4,5] 

Theorem  1 

Let  the  discrete-time  process  t Xn J n > -j  be  regenerative  with  respect 
to  the  renewal  process  Also,  let  =  Ri+i  -  R^,  i=l,2,..., 

denote  the  length  of  the  i-th  regeneration  cycle,  and  let  f  be  a  nonnegative, 
real  valued,  measurable  function. 

If  C  =  <  »  and  S  =  E{  I  f(X±)j  <  “,  then, 

i=l 

1  n  1  ,  n  S 

X.im  -  Z  f(X^)  =  £im  -  E  \  Z  f(X^)  )  =  -  ,  w.p.  1 
n+o0  n  i=1  n-«°  n  i=i  C 

Furthermore,  if,  in  addition  to  the  finiteness  of  C  and  S,  the  distribution  of 
q  is  not  periodic,  then  X^  converges  in  distribution  to  a  random  variable 
Xa,  and 

S 

E{f(X»);  =  - 
C 


Thus,  under  the  conditions  stated  above,  the  limiting  (expected)  average, 
and  the  mean  of  the  limiting  distribution  of  lf(Xn))n>^  exist,  coincide, 
and  are  finite.  Moreover,  their  common  value  is  then  given  in  terms  of  the  per 
cycLe  quantities  S  and  C. 
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Given  a  RMAA,  let  be  the  process  of  interest  associated  with 

the  random-access  system;  this  process  might,  for  example,  be  the  delay  process 
induced  by  the  algorithm.  Then,  provided  that  {xn)n>i  can  be  shown  to  be 
regenerative,  the  regeneration  theorem  itself  shows  the  way  to  establish  the 
existence  of  steady  state,  and  to  compute  the  steady-state  moments,  and  the  distribu¬ 
tion  of  { } n  ^.j,  by  appropriately  selecting  the  function  f. 

In  virtually  all  existing  RMAAs,  it  is  relatively  easy  to  identify  regenera¬ 
tive  times  (e.g.,  when  the  system  becomes  empty,  or  when  an  appropriate  Markov 
chain  hits  a  suitable  fixed  state),  at  which  the  process  of  interest  probabilis¬ 
tically  restarts  itself.  Given  a  RMAA  and  a  function  f,  the  problem  then  is  to 
exploit  the  dynamics  of  the  algorithm,  to  find  those  per  cycle  properties  of  the 
sample  function  of  the  process,  that  could  be  subsequently  used  to  evaluate  the 
quantities  0  and  S. 

In  section  3,  it  is  shown  that  for  the  delay  process,  and  for  f(x)  =  x, 
the  computation  of  S  and  C  are  intimately  related  to  the  solution  of  an  infinite 
dimensional  system  of  linear  equations.  It  can  be  shown  that  this  is  the  case 
when  f(x)  =  xn,  n  =  2,3,...,  as  well  [12].  Therefore,  the  steady-state  moments 
of  the  delay  process  induced  by  a  particular  algorithm,  can  be  computed 
from  the  solution  of  the  corresponding  infinite  linear  system.  In  Appendix  A, 
we  give  a  number  of  general  results,  that  are  useful  in  establishing  the 
existence  and  uniqueness  of  a  solution,  and  in  developing  approximations  to 
the  solution  of  such  systems.  In  section  4 ,  we  apply  these  results  to  the 
specific  infinite  linear  systems  developed  for  the  three  algorithms  of  section 
3.  This  procedure  involves  the  following  steps. 

Step  1  Find  conditions  under  which  the  infinite  linear  system  has  a 
unique,  nonnegative  solution. 


.*>  •-»*>  '-*  t-  ••  -  -■  --  -•  ■  r..  . . 

b 

Step  2  Show  that  the  variables  of  interest  coincide  with  the  unique 
solution. 

Step  3  Develop  arbitrarily  tight  upper  and  lower  bounds  on  the  solution. 

3.  THREE  ALGORITHMS  AND  THEIR  RELATED  SYSTEMS  OF  EQUATIONS 

For  all  three  algorithms  of  this  section,  we  assume  a  collision-type,  packet 
switched,  slotted,  broadcast  channel.  The  channel  is  accessed  by  a  very  large 
(effectively  infinite)  number  of  identical,  independent,  packet-transmitting, 
bursty  users.  The  cumulative  packet  generation  process  is  modelled  as  a  Poisson 
process,  with  intensity  A  packets  per  slot.  However,  the  proposed  method  can  be 
applied  equally  well,  when  the  number  of  packets  per  slot  are  independent  and 
identically  distributed  (i.i.d)  random  variables. 

We  define  the  delay,  Vn,  experienced  by  the  n-th  arrived  packet,  as  the  time 
difference  between  its  arrival  at  the  transmitter,  and  the  end  of  its  successful 
transmission.  We  are  interested  in  evaluating  the  steady  state  statistics  of 
the  delay  process  f^n^n»l'  when  they  exist.  Due  to  space  limitations  in 
this  paper,  we  give  explicit  results,  only  for  the  first  moment  of  the  delay 
process.  However,  higher  moments  of  the  delay,  as  well  as  other  quantities  of 
interest  can  be  computed,  using  the  same  method;  (the  computation  of  the  delay 
variance  for  the  third  algorithm  in  this  paper  can  be  found  in  [12]). 


3. 1  Example  1  :  Controlled  ALOHA 

The  earliest  and  most  well  Known  RMAAs  belong  to  the  class  of  the  ALOHA 
techniques  [14,6,19].  Here,  we  analyze  a  version  of  the  slotted  ALOHA 
algorithm,  that  operates  with  each  user  transmitting  a  newly  arrived  packet, 
in  the  first  slot  after  its  arrival.  Should  this  cause  a  collision,  each 


involved  user  independently  retransmits  its  packet  in  the  next  slot, 
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with  probability  f. 

A  packet  whose  transmission  is  unsuccessful  is  said  to  be  blocked.  Let  be 
the  number  of  blocked  packets  at  the  beginning  of  slot  i  (time  segment  (i,i+l] ). 
This  number  will  be  referred  to  as  the  backlog  size.  Also,  let  R^  denote  the 
number  of  blocked  packets  retransmitted  in  slot  i,  and  denote  the  number  of 
new  packets  transmitted  in  slot  i.  Given  =  m,  then  clearly, 

A 

P(Ri=r)  =  B^(f)  =  (™)  fr  (l-f)m_r  ,  i  =  0,1,2,...  (1) 

A  e_A  ln 

P(N^=n)  =  p^j  —  ——————  ,  i  —  0,1,2,...  (2) 

n! 


The  delay  process  induced  by  the  above  algorithm  "probabilistically 
restarts  itself"  at  the  beginning  of  each  slot  Tif  at  which  MT  =  0,  i  =  1,2,...; 
this  is  so  because  the  number  of  arrivals  per  slot  is  an  i.i.d.  sequence  of 


random  variables.  Precisely,  let  T.j  =  1 ,  and  define  Ti+1  as  the  first  slot 


after  T^  at  which 


M_,  =0.  The  interval  (T4  ,  T4  .  ) ,  i  =  1,2,...,  will  be 

M+1  1  11 


referred  to  as  the  ith  session. 

Let  R^ ,  i  =  1,2,...,  denote  the  number  of  packets  successfully  transmitted 
in  the  interval  (0,  T^+^]  (Note  that  R ^  also  represents  the  number  of  packets 
arrived  during  the  interval  (0,  T^-1)).  Then,  C  ^  -  R i  =  1,2,...,  is 

the  number  of  packets  successfully  transmitted  in  the  interval  (TL,  Tl+^;.  I’he 
sequence  is  a  renewal  process,  since  is  a  sequence 

of  nonnegative  i.i.d.  random  variables.  Furthermore,  the  delay  process 
is  regenerative  with  respect  to  the  renewal  process  (R  i  ,•) ,  with  regenera¬ 
tion  cycle,  C i . 

Ftom  theorem  1,  with  f(D^)  -  V  ,  we  have  that  it  C  -  Kit""  ^ and  it 


vvv 


£ {  £1?^)  <  ■*>,  then,  there  exists  a  real  number  D,  such  that, 
i=1 


1  n  In  3 

D  =  Him  -  ZV^  =  Jtim  -  E{  ZV^ )  =  E{  P  »)  =  -  a.e. 

n-Hx,  n  t=i  n-*°  n  i=1  C 


Thus,  provided  that  both  S  and  C  are  finite,  the  limiting  average,  the 
limiting  expected  average,  and  the  mean  of  the  limiting  distribution  of  i > 1 » 

exist,  coincide,  and  are  finite;  their  common  value  D  will  be  referred  to  as  the 
mean  packet  delay. 

Next,  we  develop  two  systems  of  equations,  whose  solution  .ay  be  used  to 
compute  the  mean  cycle  length  C,  and  the  mean  cumulative  delay  S.  The 
properties  and  the  computation  of  the  solution  will  be  postponed  until  section  4. 


l.a  Mean  Cycle  Length 

If  the  mean  session  length  H  =  E  (Ti+i  -  T^i  ,  i  >  1,  is  finite,  then 
by  Wald's  identity,  we  have  that. 


c  =  ah 


To  determine  H,  we  proceed  as  follows.  bet  hx  denote  the  random  number  of 
slots  needed  to  return  to  zero  backlog  size,  starting  from  a  slot  j  where  the 
backlog  size  is  equal  to  i,  i  ?  0.  By  definition,  then,  4-3  ihQJ.  The 
operation  of  the  algorithm  yields  the  following  relation  for  the  h^'s. 


ho  \ 


1  +  hK 


if  =  0,  1 


if  Nj  >  1 


(4. a) 


1  +  h. 


if 

*3 

N3 

=  0 

if 

R. 

+ 

N  ■ 

=  1 

3 

3 

if 

R-i 

N_i 

>  1 

3 

3 

(4 .  b ) 


If  we  let  =  E  i hx ] ,  1  <>  0 ,  then  after  taking  expectations  in  (4)  we  obtain, 


Hi  =  bi  +  Z  cik  Hk 


where  bj.  =  1,  i  >  0,  c0o  =  ^01  =  b,  c0i  =  py ,  i  *  2,  cik  =  Pk-i*  k  >  i  +  1, 

Ci,  i+1  =  p1  (1-B0(f))'  Cii  =  P0(1-B1(f))  +  ?1B0(f)'  ci,i-1  =  P0B1(f)' 

=  0,  k  <  i-1 ,  and  where  p^,  B^(f),  i  >  0,  0  <  j  4  i,  are  as  defined  in 

(1),  and  (2),  respectively. 

Note  that  the  mean  session  length  H,  can  be  computed  from  system  (5),  since 
H  =  H0. 

l.b  Mean  Cumulative  Delay 

The  mean  cumulative  delay.  S',  can  be  computed  using  a  system  of  equations 
similar  to  system  (5).  To  develop  such  a  system  we  proceed  as  follows.  Let  Wy 
denote  the  cumulative  delay  experienced  by  all  the  packets  that  were  successfully 
transmitted  during  hx  slots.  also,  let  =  Elw^},  i  *  0,  and  note  that 

c  =  w0. 

The  operation  of  the  algorithm  yields  the  following  relations  for  the  w^'s. 
N-j  if  Nj  =  0,  or  1 

wo  {  <°-a> 


if  N j  =  0 ,  or  1 


(O.a) 


H3  +  WN. 


if  N-j  >  1 


(1)  Here,  for  convenience,  we  count  the  delay  of  a  packet,  starting  from 
the  bejinning  ot  the  first  slot  after  its  arrival. 
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( 1 

+ 

+ 

W  • 

1 

if 

*3 

+ 

N3 

=  0 

H 

II 

+ 

Nj 

+ 

Wi+N  ,-1 

if 

rd 

+ 

N3 

=  1 

(6.b) 

i>0  1 

3 

\  i 

Nj 

Wi+N  . 

1 

if 

*3 

+ 

nd 

>  1 

After  taking  expectations  in  (6)  we  obtain 


wi 


*  bi  + 


'ik 


W,. 


k=0 


(7) 


where  b{j  =  ^  ,  b^  =  i  +  a  ,  i>  1,  and  ci)c  are  as  defined  in  (5). 

3.2  Example  2  :  The  "0.487"  Algorithm 

This  algorithm  is  the  most  efficient  RMAA  known  to  date,  for  the  Poisson 
infinite-user  population  model  and  ternary  feedback;  (it  allows  a  maximum 
throughput  of  0.487  packets  per  slot).  It  is  assumed  that  at  the  end  of 
each  slot  i,  the  users  receive  a  feedback  =  0,1,  or  c,  if  in  slot  i  there 
were  zero,  one,  or  more  than  one  packets  transmitted. 

The  following  is  a  brief  description  of  the  algorithm;  for  details, 
motivation,  and  background  discussions,  the  reader  is  referred  to  [8],  and 
[18]  . 

Suppose  that  at  the  beginning  of  slot  v(time  segment  (v,v+l]),  all  packets 
that  arrived  oefore  time  tv  <  v,  have  been  successfully  transmitted,  and  there 
is  no  information  concerning  the  packets  that  may  have  arrived  in  the  interval 
[ ty , v ) ,  (i.e.,  the  distribution  of  the  interarrival  times  of  the  packets  in 
[tv,v)  is  the  same  as  the  one  assumed  originally).  The  beginning  of  such  a 
slot  v  is  calLed  a  "collision  resolution  instant".  The  time  difference 
dv  -  v  -  tv  wiLl  be  referred  to  as  the  "lag  at  v".  In  slot  v,  the  users 
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that  generated  packets  in  the  interval  [tv,  ty  +  Uv),  where  Uv  =  min(dy,  A), 
are  allowed  to  transmit;  A  is  a  parameter  to  be  properly  chosen  for  throughput 
maximization.  In  this  case,  we  say  that  the  interval  [tv,  tv  +  Uv)  is  "trans¬ 
mitted".  Depending  on  the  received  feedback  zv,  the  algorithm  operates  as 
follows: 

If  zv  =  0,  or  1,  the  transmitted  interval  is  "resolved",  and  at  time 
v  +  1  the  resolution  of  another  interval  starts,  where  now  tv+1  =  tv  +  Uy; 
v  +  1  is  a  new  collision  resolution  instant. 

If  zv  =  c,  the  collision  in  the  transmitted  interval  is  resolved  according 
to  the  following  steps: 

1.  The  interval  [ty,  ty  +  Uy/2)  is  transmitted  in  slot  v  +  1. 

1.1  If  zv+1  ~  c,  the  resolution  of  the  interval  [ty,  tv  +  Uv/2)  starts 
at  v  +  2  according  to  step  1,  with  tv+2  =  tv,  and  Uv+2  =  Uv/2.  The 
users  in  [tv  +  Uv/2,  ty  +  Uy)  stop  participating  in  the  collision 
resolution  process  (in  this  case  we  say  that  the  interval 

[tv  +  Uy/2,  tv  +  Uy )  "returns"  to  the  unexamined  portion  of  the 
arrival  axis). 

1.2  If  Zy+i  =  0,  the  resolution  of  the  interval  [tv  +  Uv/2,  tv  +  Uv) 
starts  at  v  +  2  according  to  step  1,  with  ty+2  =  ty  +  Uv/2,  and 

Uv+2  =  V2' 

1.3  If  zv+i  =  1,  the  interval  [tv  +  Uv/2,  ty  +  Uy)  is  transmitted  in 
slot  v  +  2. 

1.3.1  If  zv+2  =  c,  the  resolution  of  the  interval  [tv  +  Uv/2,  tv  +  Uv) 
starts  at  v  +  3  according  to  step  1,  with  tv+j  =  ty  +  Uy/2,  and 

Uv+3  =  uv/2* 


1 1 

1.3.2  If  zv+2  =  the  originally  transmitted  interval  has  been  resolved, 
and  the  resolution  of  another  interval  starts  at  v  +  2,  with 
tv+2  =  tv  +  Uv/2;  v  +  2  is  a  new  collision  resolution  instant. 

Et>r  the  analysis  of  the  above  algorithm,  we  need  the  following  definitions 

6  :  the  length  of  the  originally  transmitted  interval  that  is  not 

returned  to  the  unexamined  portion  of  the  arrival  axis;  we 
refer  to  6  as  the  "examined  portion  of  Uv" 
i  :  the  number  of  slots  needed  for  the  resolution  of  the  interval 

[tv,  tv  +  &). 

N  :  number  of  packets  in  6 

u)  :  sum  of  delays  of  the  N  packets,  after  the  resolution  process 

of  [tv,  tv  +  Uv)  begins 

y  :  sum  of  delays  of  the  N  packets,  until  the  instant  ty  +  Uy. 

E{x|u}  :  conditional  expectation  of  the  random  variable  X,  given  that 

Uv  =  u 

Let  be  the  sequence  of  successive  collision  resolution 

instants, and  let  d^  be  the  lag  at  v^.  It  is  known,  [10],  that  the  sequence 

1 ' ■* i ;  i  *  1  LS  a  Mar*<ov  chain,  with  state  space,  F  a  denumerable  dense  subset 

of  the  interva]  [1,  *) .  Let  -  1,  =  1,  and  define  T  +^,  as  the  first  slot 

after  Titat  which  dT  =  1.  From  the  description  of  the  algorithm  it  can  be 

i+l 

seen,  after  a  little  thought,  that  the  induced  delay  process  probabilistically 
restarts  itself  at  the  beginning  of  each  slot  T^,  i  =  1,2,...  .  Therefore, 
using  the  notation  and  definitions  of  example  1,  the  mean  packet  delay  D  is 
equal  to  .-l/C  provided  that  both  H  and  C  are  finite. 
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2 .a  Mean  Cycle  Length 

As  in  example  1,  if  the  mean  session  length  H  =  E{'P^+1  -  T^l  is 
finite,  then  C  =  Ah.  To  evaluate  H  we  proceed  as  follows. 

Let  h<j  denote  the  random  number  of  slots  needed  to  return  to  lag  equal 
to  one,  starting  from  a  collision  resolution  instant  with  d^  =  d.  Note  that, 
by  definition,  h^  is  the  session  length.  The  operation  of  the  algorithm  yields 
the  following  relations  for  the  hd's,  d  t  F. 


1  <,  d  < 


A  ,  hd  =  | 


i  if  Z  =  1 

fc+hd-6+2  if  ^  >  1 


(8. a) 


d  >  A  ,  hd  =  «-+hd„5+fc 


(3.b) 


Taking  expectations  in  (8)  yields: 


1  <  d  <  A  ,  Hd  =  El£|d}+Z  p(r,s|d)Hd_r4.s 

r,s 

s;*1 


(9. a) 


d  >  A  ,  Hd  =  E l £ | A}  +  L  p(r,s | A)Hd_r+s 

s,r 


(9.b) 


where  p(r,s|x)  is  the  joint  conditional  probability  distribution  of  6,  and 
i,  at  the  point  values  r  and  s,  given  that  the  transmitted  interval  is  of 
length  x.  Note  that. 


( 1+Ax)e 


■Ax 


p(r ,1 |x)  = 


if  r  =  x 
otherwise 


System  (9)  can  be  written  in  the  form 


Hd  “  bd  +  ^  °dt  Ht  »  ,lu  r 

f.  F 


(  10} 


where  bd  =  E(£|d),  1  ^  d  <  A,  bd  =  E  t  £  |  A }  ,  d  >  A,  and  where 

cdt'  c  e ^  are  nonnegative  coefficients  that  can  be  appropriately  identified 
from  (9).  The  conditional  expectation  E{£)d},  1  =*  d  <  A,  can  be  computed 
as  shown  in  Appendix  B. 

2.b  Mean  Cumulative  Delay 

Let  wd  denote  the  cumulative  delay  experienced  by  all  the  packets  that 
were  successfully  transmited  during  hd  slots.  The  operation  of  the  algorithm 
yields  the  following  relation  for  the  wd's,  d  e  f . 

!UJ+<P  if  £  =  1 

(10. a) 

w+y+wd_d+x  if  £  >  1 


iuJ  +  y+(d-A)N  if  £  =  1 

(10. b) 

u)+ip+( d-A)N+wd_,5+£  if  £  >  1 

Taking  expectations  in  (10),  yields: 

1  <  d  <  A,  Wd  =  E (u)  |  d}  +  Ei4/|d}  +  l  p( s , r | d ) Wd_r+S  (11. a) 

r  ,s 
s^1 

d  >  A,  Wd  =  eN|j/  +E 1  V  |  a)  +(  d-  i )  E{N  |  A/  +  l  p( s , r | A) Wd_rfs 

s,r 

(11. b) 


System  (11)  can  be  written  in  the  form 


l2 
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where  b£=E{w  | di+2 1  y  | d)  ,  L  >>  d  <  A,  b£=Eiw  |  n}+El  v|  d}  t  ( d-A) E l N  |  &}; 
and  where  the  coefficients  cdt,  d,  t  £  F  are  as  defined  in  (10).  The 
conditional  expectations  Eiwjdi,  E{</|d},  1  <  d  4  A,  and  E { N |  A>  can  be 
computed  as  shown  in  Appendix  B. 

3.  Example  3  :  Stack  Algorithm 

A  new  trend  towards  the  design  of  easy-to-implement  RAAs,  which  combine 
stability  and  good  performance,  with  modest  feedback  requirements,  started  with 
the  introduction  of  the  "stack"  algorithm  by  Tsybakov  and  Vvedenskaya  [9] .  The 
new  class  of  algorithms  has  "limited  feedback  sensing"  and  "free  access" 
characteristics.  Limited  feedback  sensing  algorithms  require  that  users  sense 
the  feedback  broadcast  only  while  they  have  a  packet  to  transmit,  and,  there¬ 
fore,  they  have  practical  advantages  over  continuous  feedback  sensing  algor ithms, 
such  as  the  algorithm  in  example  2.  The  "free  access"  characteristics  of  the 
new  algorithms  simplify  their  implementation,  since  newly  arrived  packets  are 
transmitted  in  the  first  slot  after  their  arrival. 

The  applicability  of  the  proposed  method  to  the  analysis  of  algorithms 
from  this  new  class  has  been  demonstrated  in  [12] ,  where  a  representative 
algorithm,  called  n-ary  stack  algoritrhm  ( SAn ) ,  has  been  analyzed.  The 
analysis  presented  in  [ 1 2 j  is  included  here  for  reasons  of  completeness. 

The  3An  uses  binary  feedback  of  the  "collision-No  Collision"  ( C-NC)  type; 
that  is,  at  the  end  of  each  slot,  the  users  that  monitor  the  channel  are 
informed  whether  that  slot  contained  a  collision  or  not.  Let  z^  denote  the 
binary  feedback  corresponding  to  slot  k,  where  zk  =  NC,  and  zk  =  C  represent 
respectively  a  noncollision  versus  collision  slot,  k.  Let  some  packet  arrive 
during  slot  k-1.  We  then  assume  that  the  packet  attempts  transmission  during 
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slot  k,  and  it  observes  the  feedback,  zk  and  all  the  feedbacks  after  that, 
until  successfully  trnasraitted.  We  also  assume  that  the  packet  has  no  know¬ 
ledge  of  the  channel  feedback  history,  z^:  i<k  (limited  feedback  sensing).  In 
its  effort  to  be  transmitted  successfully,  the  packet  utilizes  a  counter,  whose 
indication  at  the  beginning  of  slot  k  is  denoted  by  1^,  and  it  applies  the 
following  set  of  rules. 

1.  A  packet  arrived  during  slot  (k-1),  sets  1^  =  1 

2.  A  packet  attempts  transmission  within  slot  k,  if  and  only  if,  Ik  *  1. 

The  packet  is  successfully  transmitted  within  slot  k,  if  and  only  if, 

1^=1  and  zJc=NC. 

3.  The  updating  of  the  counter  indication  Ik  is  done  as  follows. 

a)  If  zk=NC,  and  IkJ*l,  then,  I^^I^-l 

b)  If  zk=C,  and  1^=1,  then,  I^+i^Jf  where  J  is  an  integer  random  variable, 
uniformly  distributed  on  (l,2,...,n),  and  n  is  an  integer  parameter, 

n  ^  2. 

c)  If  Zk=C,  and  I)t>l,  then  I^+^I^t  -1 

The  operational  characteristics  of  the  SAn  are  perhaps  better  explained  by 
introducing  the  concept  of  a  "stack",  as  it  was  first  done  in  [9).  A  stack  is 
an  abstract  storage  device,  consisting  of  an  infinite  number  of  cells,  labelled 
1,2,3...  .  The  number  of  packets  that  a  cell  can  accommodate  is  unrestricted. 

At  the  beginning  of  each  slot,  k,  the  m-th  cell  of  the  stack  contains  the 
packet  with  I)c=m.  Packets  are  eventually  successfully  transmitted  after  moving 
through  the  cells  of  the  stack  in  accordance  with  the  rules  defined  above,  as 
shown  in  figure  1.  To  resolve  conflicts,  the  algorithm  splits  uniformly  the 
group  of  collided  packets  into  the  first  cells  of  the  stack.  The  integer  n 
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is  an  algorithmic  parameter,  whose  value  may  be  chosen  for  performance 
optimization. 

The  random  access  system  operates  with  the  SAn  in  sessions.  A  session  is 
a  sequence  of  consecutive  slots  that  begins  and  ends  at  two  consecutive 
algorithmic  renewal  instants,  to  be  defined  below.  Those  instants  are  denoted 
by  T^,  i  >  1,  and  are  determined  by  means  of  a  renewal  counter.  The  first 
session  begins  at  T^=1,  with  the  renewal  counter  set  to  "2".  Depending  on  the 
channel  outcomes,  the  renewal  counter  updates  its  indication  in  accordance  with 
rules  3. a,  and  3.c.  The  second  renewal  instant,  T2,  is  the  instant  at  which 
the  renewal  counter  drops  to  “1"  for  the  first  time;  this  signifies  the  end  of 
the  first  session.  Immediately  after  T2  the  renewal  counter  is  reset  to  "2", 
and  the  second  session  begins.  This  process  continues  indefinitely,  and  defines 
the  successive  sessions  [Tj_ ,  Ti+j),  i=1,2,...  . 

From  the  definition  of  the  session  given  above,  it  can  be  easily  seen  that 
immediately  before  the  end  of  a  session,  all  cells  of  the  stack  are  empty  of 
packets,  and  that  a  new  session  begins  with  the  group  of  new  packets,  that 
arrived  during  the  last  slot  of  the  previous  session,  placed  in  cell  #1.  Due  to 
the  independent  and  stationary  increments  property  of  the  arrival  process  the 
session  lengths  T^+i-Ti,  i=1,2,...,  are  positive  i.i.d.  random  variables. 
Furthermore,  as  in  the  previous  two  examples,  the  delay  process 
is  regenerative  with  respect  to  the  renewal  process  ,  with 

regeneration  cycle  length,  the  number  of  packets  successfully  transmitted  during 
a  session.  Therefore,  the  mean  packet  delay  D  induced  by  the  SAn  is  equal  to  S/C, 
provided  that  both  S  and  C  are  finite. 


Mean  Cycle  Length 


As  in  the  previous  examples,  c  =  ah. 


To  determine  the  mean  session 


length,  H,  we  proceed  as  follows.  Let  h^,  i  >  0,  denote  the  random  length  of 
a  session  that  begins  with  i  packets  in  the  first  cell  of  the  stack.  After 


a  little  thought,  it  can  be  seen  that  the  dynamics  of  the  SAn  yield  the 
following  relation  for  the  h^'s: 


n 

hg  =  h i  =  1  ;  hi  =  1  +  E  hj  ,  i  >  2  (13) 

j=l  3  3 


where  N^,  N2»«.*,Nn  are  independent,  Poisson  random  variables  with  parameter 
A,  which  are  also  independent  of  the  random  variables  I1#  I2,  I...,  In, 
which  are  multinomially  distributed,  with 


Pdfi 


,]tn  ^-n^ 


il 

i  1  * • • • in I 


n 

0  <  ij  <  i,  2  ij 
j=1 


i 


A 

If  we  define  =  Eth^}  =  E{hI|l=i},  then,  after  taking 
expectations  in  (13),  we  obtain. 


Hg  =  H,  =  1 

00 

“i  =  bi  +  1  cik  Hk  ,  i  >  2  (14) 

k=2 

where  ^  =  1  +  c±1  +  ciQ,  i  >  2,  cik  =  n(pk*B^(n-1 ) ) ,  i  >  0,  k  >  0 
and  where  B^C),  pk  are  as  defined  in  (1),  and  (2),  respectively,  and 
*  signifies  convolution. 

The  mean  session  length  is  given  by 
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3.2  Mean  Cumulative  Delay 

Let  w^  ;  i  >0,  denote  the  cumulative  delay  experienced  by  all  the 
packets  successfully  transmitted  during  a  session  that  begins  with  i  packets 
in  the  first  cell  of  the  stack.  The  rules  of  the  algorithm  yield  the  following 
relation  for  the  w^'s: 

n- 1  n 

wQ  =  0;  w1  =  1;  wi  =  i  +  £  Q.  hj  +N  +  £  Wj  +N  ,  i>1  (16) 

3  =  1  j  3  3=1  3  3 

where  I j,  N j ,  i  <  j  <  n.  are  as  defined  in  (13)/  and 

n 

Qj  -  E  I*'  1  <  3  <  n-1 
m=  j+1 


Taking  expectations  in  (16)  yields  the  following  system  of  equations  for  the 

A 

conditional  mean  cumulative  delay  =  Etw^}  =  E(wI|l=i},  i  >  0: 


Wn  =  0,  W.  =  1 


wi  =  bi  +  Z  cik  wk'  1  '  2 
k=2 


(17) 


where 


n-2  i  00  i-j  .  /  £  \  1  \ 

bi=bi( {Hi}i  *))=i+ci1+  4  L  L  1  Bj  '  ®k  ( - Pm(i-3-l<)r 

4=0  ]  =  0  m=0  k=0  '  n  ’ 


rtk+m' 1  '  2' 


and  where  c^,  b£(  •  )  ,  p^ ,  i  *2,  k  *2  are  as  defined  in  (14),  (1),  and 
(2),  respectively. 

The  mean  cumulative  delay,  S,  is  given  by 


oo 


1=0 


s 


( lb) 
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4 .  SYSTEM  SOLUTION  AND  MEAN  PACXfiT  DELAY  BOUNDS 

In  this  section,  we  investigate  the  conditions  under  which  the  infinite 
dimensional  linear  systems  (5),  (7),  (10),  (12),  (14),  and  (17)  have  unique, 
nonnegative  solutions,  and  we  develop  upper  and  lower  bounds  on  those  solutions. 
These  bounds  are  then  used  to  obtain  bounds  on  the  mean  packet  delay.  We  pro¬ 
ceed,  following  the  steps  outlined  in  section  2. 

4 . 1  Step  1 

Etor  convenience,  we  rewrite  an  infinite  linear  system  in  an  operator  form. 
Specifically,  let  E  be  the  space  of  sequences  X  =  {x(v)i:  A  ♦R*  where  A 
is  a  countable  set.  Also,  let  E^  be  the  subspace  of  E  for  which, 

£  |cMvx(v)  |  <  «,  utA 
veA 

We  define  the  operator  L  =  (Lp(x)}:  E^-^E,  as  follows. 

Vx)  =  \L  +  Z„cLuvx(v)  '  *eA  '  x  e  fL 

VE A 

In  this  notation,  systems  (5),  (7),  (10),  (12),  (14),  and  (17),  can  be  written 
in  the  form, 

SL  =  L(SL),  Sl  e  El  (19) 

We  are  interested  in  the  existence  and  uniqueness  of  nonnegative  points 
S*"1  *■  EL,  that  satisfy  (19);  such  points  will  be  referred  to  as  fixed 
points  of  L,  and  represent  solutions  to  the  corresponding  infinite  linear 
system  of  equations.  The  question  of  uniqueness  of  a  fixed  point  SL,  or 
equivalently  of  the  solution,  iS^(i)},  to  the  system  that  operator  L 
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represents,  depends  upon  what  conditions  are  imposed  on  the  solution.  Thus, 
after  the  existence  of  a  solution,  (s^Ci)},  has  been  established,  one 
has  to  indicate  a  class  of  sequences  in  which  the  solution  is  unique.  If  the 
algorithmic  sequences  of  interest  or  {w^l  belong  to  the  indicated 

class,  then  they  must  coincide  with  the  solution  is^(i)i.  (This  will  be 
examined  in  Step  2 ) . 

Appendix  A  includes  a  number  of  results  that  can  be  used  to  establish 
existence  and  uniqueness  of  a  fixed  point  of  an  operator.  Depending  on  the 
operator,  some  are  more  straightforward  to  apply  than  others.  Among  the 
results  in  Appendix  A  that  can  be  used  to  establish  existence  of  a  solution. 
Lemma  A. 2  is  usually  the  most  useful.  According  to  Lemma  A. 2,  to  establish 
existence  of  a  nonnegative  fixed  point,  SL,  it  suffices  to  find  a  point 
x°  e  El,  such  that, 

0  <*  L(X°)  <  X°  (20) 

A  point  X°,  satisfying  (20),  also  serves  as  an  upper  bound  on  SL.  Further¬ 
more,  to  establish  a  lower  bound  on  SL,  it  suffices  to  find  a  point  Y°  t 
such  that, 

Y°  *  L(Y°)  <■  X°  (21) 

Thus,  under  (20)  and  (21),  we  have  that, 

Y°  *  SL  4  x°  (22) 


We  proceed  now  with  the  analysis  of  the  systems  developed  in  section  3. 


21 


1.  Controield  ALOHA 

System  (5)  —  Existence;  System  (5)  corresponds  to  an  operator  Li  with 
Li  LL 

bU  =  ,  cMV  =  CjIV ,  u ,  v  e  Nq,  where  NQ  is  the  set  of  nonnegative  integers, 

and  the  b^’s  and  c^v's  are  as  defined  in  (5).  If  we  let  X°  =  ix°(k)i  with 
x°(k)  =  ot^k  +  p^,  k  A  o,  then  by  straightforward  manipulations  we  have  that,  for 
this  choice  of  XQ,  (20)  is  satisfied  if  and  only  if  the  following  inequalities  are 
satisfied. 

9  A  <  4k(f)  -  pQB^(f)  +  pj^B^ff),  for  every  k  A  1  (23) 

«u  A  Sup| -  ,  k  A  1  >  (24  ) 

Uk(f)-A  J 

|JU  *  (1  +  VA  -  Pl>>/<Po  +  Pi)  (25) 

It  can  be  readily  seen  from  (23)  that  if  the  retransmission  probability  f  is 
constant  in  every  slot,  then  there  is  no  A  A  0  for  which  (23)  is  satisfied. 

If  the  retransmission  probability  f^,  at  each  slot  i,  were  allowed  to  depend 
on  the  current  backlog  size,  ,  in  accordance  to  a  stationary  control  policy 
f  =  f(f“r),  then  it  is  of  interest  to  choose  f(  )  so  that  it  maximizes  the  set 
of  A's  for  which  inequality  (23)  is  satisfied.  This  is  equivalent  to 
maximizing  .,k(f)  with  respect  to  f .  It  can  be  easily  verified  that,  for 
every  k  A  1,  -(f)  is  maximized  for  f(k)  -  t*(k),  where  (^) 

2.  We  should  mention  that,  in  a  distributed  environment,  the  backlog  size 
dependent  retransmission  probability  f*(*)  is  nonimplementable ,  since 
users  are  not  aware  of  the  current  backlog  size.  However,  the  control 
policy  given  by  (26)  can  be  implemented  approximately  by  adaptive  control 
schemes  that  estimate  the  current  backlog  size  using  observable  feedback 
information  from  the  past  activity  on  the  channel  [6,  19,  20]. 
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1-A 

c*(k)  =  —  ,  k  »  i 
k-A 


From  this  point  on,  we  assume  that  f  is  chosen  as  in  (26).  Under  this  assumption. 


inequality  (23)  is  satisfied,  provided  that. 


a  <  inf  Uk<£*) ,  k  i  1}  =  e_1 


To  satisfy  inequalities  (24),  and  (25),  we  choose. 


Ou  =  sup 


ok(f*)-A 


•  '!•  s 


l+iu(A-Pl)  1 

'  ^U  =  T"7 

A  e_  ~A 


Similarly,  it  is  straightforward  to  show  that  if  V  <  e-1  then  the  point  Yc 
with  y°(k)  =  a^k  +  1^,  k  >0,  and 


H  =  ^ 

e_1-A 


1-A  1 

*+*  I'"-  A 


satifies  (21).  Thus,  from  (22)  and  for  A  <  e“  we  have  that  system  (5)  has 

Li  L  i 

a  solution,  S  =  Is  (k)},  such  that 


0  <  u^k  +  *  s  ( k )  <  uuk  +  ou  ,  k  *  0 


where  uu,  6U  are  as  given  by  (27),  and  i K ,  are  as  given  by  (28). 


System  (7)  —  Existence:  System  (7)  corresponds  to  an  operator  L2  with  bp  =  b^ * , 

L2 

C|JV  =  c[lv,  ,  v  t  N0,  where  the  b^'s  and  cMV's  are  is  defined  in  (7). 


Due  to  the  fact  that  bk  is  a  linear  function  of  k,  and  since 


1,  i 


1 
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it  can  be  easily  seen  that  there  is  no  linear  sequence  X°  =  {x°(k)i 
satisfying  (20).  However,  given  A  <  e-^ ,  it  is  straightforward  to  show  that 
we  can  choose  coefficients  Yu»  6U,  £u,  Y£»  <$£»  4£  such  that  the  point  X°  with 
x°(k)  *  Yuk2  +  <$uk  +  Su»  k  >  0,  and  the  point  Y°  with  y°(k)  =  6yk2  +  6£  + 
k  >  0,  satisfy  (20),  and  (21),  respectively.  The  following  is  such  a  choice: 

.  A+y  (A2+A+e~1)  1+Y  (l+A-e“A)+6  (1-e-A) 

1  x  u  A  u  u 

Yu  —  ,  ou  —  ,  Cu  ~  A  " 

2(e“1-A)  e_1-A  ( 1+A)e'A 


1  A+YA(A2+A+e"A(l-2A)) 

Yg,  =  -  /  - - - - 

2(e-A-A)  e-A-A 


H  =  x 


l+Y£(l+A-e-A)+6^(l-e“A) 

(l+A)e“A 


h2  1-2 

Then,  if  A  <  e”1,  then  system  (7)  has  a  solution  S  =  {s  (k)}  ,  such  that. 


0  <  Yskz  +  6jk  +  Cj <  s  (k)  <  6uk^  +  6uk  +  Cu  #  k  >0  (32) 


where  Yu,  <SU,  Cu  are  as  given  by  (30),  and  Yjj,,  <$£,  <»£  are  as  given  by  (31), 


Systems  (5)  and  (7)  —  Uniqueness 

L1 

We  will  show  that  both  the  solution  (s  (i)i  of  system  (5)  and  the 
L2 

solution  is  (i)l  of  system  (7)  are  unique  in  the  class 


h 2  =  \  X  :  sup 


|x(i) 


itNQ  i^+c 


where  c  is  a  positive  constant. 

We  start  with  system  (7).  Since  L2  is  majorant  of  itself,  from  theorem  A. I, 


we  have  that  L2  has  a  principal  fixed  point  S*  ,  such  that  0  *  S*  <  S 

l2 

According  to  theorem  A. 2,  the  fixed  point  S  is  unique  in  the  class, 


I*2  (  |x(i)  |  ) 

I  sL2  ) 


^2  l2 
provided  that  Y °  £  E*  .  Since,  by  definition,  r  £  E2,  S  will  be  unique 

L2 

in  E2,  if  we  show  that  E*  =  E2.  According  to  lemma  A.1,  it  suffices  to  show 
that. 


s*  (i) 


i£N_  i^+c 


s*  ( i) 

inf  - >  0 

i  £N  i  2+c 

o 


u2  u2 

Since  0  <  s*  (i)  'is  (i),  (33)  follows  from  (32).  To  show  that  (34) 

L2 

holds,  we  use  the  power  sequence,  tSn  }nJ1,  of  L2  with  initial  point  0.  By 

L2 

definition  (see  Appendix  A),  Sn  is  the  point  that  results  after  L2  operates  n 

l2  l2  l2 

times  on  the  zero  point,  (i.e.,  Sn  =  L^fO)),  and  Sn  +  S*  ,  as  n  +  »  .  Due 
L2  l2 

to  the  fact  that  >0,  c^^  '0,  i,  k  £  NQ,  we  have  that  0  <  sn  (i)  *> 

l2  l2 

sn+i(i)  <  s*  (i),  for  every  n  *  1,  i  >0.  Also,  it  can  be  readily  shown  by 
induction  that,  for  every  i  J  1,  n  >  1, 


s„  ( i )  >  m  -  n( n-1 )/2 


From  ( 35 )  we  obtain 


l2  l2  2 
s*  (i)  (i)  1  i  +i  1 

£im  inf  — -  >  iim  inf -  >  -  - -  -  (36) 

i-H»  i2+c  i  -*»  i2+c  2  i2+c  2 

l2 

(34)  follows  from  (36),  and  the  fact  that  s*  (i)  >  0,  i  *0. 

L1 

The  uniqueness  in  Ej  of  the  solution  {s  (i)}  of  system  (5)  follows 
from  theorem  A. 4,  part  (ii),  after  one  identifies  L1  with  02  and  L2  with  0y, 
in  the  theorem. 

2.  The  "0.487"  Algorithm 

System  (10)  —  Existence  and  Initial  Bounds 

L1  Li 

System  (10)  corresponds  to  an  operator  with  bp  =  ,  c^v  =  Cjjv, 

JJ,  v  e  F,  where  the  b^'s  and  cMv's  are  as  defined  in  (10).  To 
establish  the  existence  of  a  nonnegative  solution  to  system  (10),  we  follow 
the  same  procedure  as  in  system  ( 5 ) . 

Let  X°  =  (x°(d) }  with  x(d)  =  otud  +  (5U,  d  eF  ,  and  let  X’  =  L^X  °). 

After  straightforward  manipulations,  we  obtain, 

x*  (d)=x°(d)+E{£|d}+au(EU|d}-E{<5|d}-(l+Ad)e-Ad-tfu(l+Ad)e-^d  ,  1  <  d  <  A 

(37. a) 

x'(d)  =  x°(d)  +  EU|d)  -  ( E  { 6 1  A}  -  eU|u}),  d  >  u  (37. b) 

According  to  Lemma  A. 2,  to  establish  the  existence  of  a  nonnegative  fixed  point 
of  L<|,  it  suffices  to  show  that  there  exist  uu,  p>u,  such  that, 


0  <  x'(d)  <  x°(d) 


for  every  d  t  F 


<3fcl) 


r.  -9:  '*1  r. ' 


TV.-  !  '  ■ 
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EYom  (37. b),  we  see  that  this  is  possible  only  if 

EU|A}  >  eU|A}  (39) 

If  (39)  holds,  then  it  can  be  readily  seen  from  (37)  that  (3a)  is  satisfied,  if 
we  choose  c^,  Pu  as  follows: 


eU|  a} 

u  E{6ja}  -  e{£.Ja} 


(40. a) 


du  =  max{-  au,  sup  (P(d))}  (40. b) 

1  <d<A 


where 

E{ A  |d>  +  au(EU|d)  -  E(6|d>  -  ( 1+Ad)exp(-Ad) ) 

p(d)  =  - 

( 1+Ad)exp(-Ad) 

The  conditional  expectations  appearing  in  the  above  expressions  can  be  computed 
as  shown  in  Appendix  B. 

Similarly,  it  can  be  shown  that,  under  (39),  the  point  Y°  =  {y°(d)> 
with  y°(d)  =  a^d  +  t d  e  F  satisfies  the  inequality  Y°  *  X,(Y°)  <  X°, 
if  and  are  chosen  as  follows: 


=  uu  '  =  in£  (  P(d)  ) 

1<d<A 


(4  1  ) 


where  «u,  P(d)  are  as  given  by  (40). 

Thus,  if  (39)  holds,  then  from  lemma  A. 2  we  have  that  system  (10)  has  a 
L1 

nonnegative  solution  S  ,  such  that. 


L1 

Ufcd  +  <  S  (d)  <  aud  +  Pu  , 


d  t  F 


(42) 
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where  u^,  tfu,  and  a^,  are  as  given  by  (40)  and  (41),  respectively. 


System  (12)  —  Existence  and  Initial  Bounds 

Let  L2  be  the  operator  that  corresponds  to  system  (12).  Also,  let 
X°  =  {x°(d)}  with  x°(d)  =  Yud2  +  6^6.  +  d  tF  ,  and  Y°  =  ly°(d)i  with 

y°( d)  =  y^d2  +  6£d  +  ,  d  e  F. 

following  the  same  procedure  as  for  system  (10),  we  can  show  that  if  (39) 

l2  l2 

holds,  then  system  (12)  has  a  nonnegative  solution  S  =  {s  (d)j,  d  £  F, 
such  that. 


Y£d2  + 


S£d  + 


^2 

<  s  (d) 


V*  +  6U<5 


+  t, 


u 


(43) 


where , 


■u 


y<L  = 


E{n| A} 


2(E{ 6 | A}  -  e{£ | A} ) 


E(w|A}  +  E { tp I  A>  -  A  e{n|A)  +  Yu  E{  (  6-i2  |  A) 
E{  6 1  A}  -  e{  i  |  A} 


Cu  =  sup  ( <M d) )  ,  (.g  =  inf  (9(d)) 
1<d<A  14d<A 


4>(d) 


2{oj  |d}+E(9|d}+Yu(E{  (5-H)2|d}-2dE{6-Ji,|d>-(l+Ad)e~XA)-6u(E{5-(i.|d}-(l+Xd)e_Xd) 

(l+Xd)e"Ad 


The  conditional  expectations  in  the  above  expressions  can  be  computed  as 


shown  in  Appendix  B. 


Remark  It  is  known  [7]  that  inequality  (39)  is  satisified  if  A  <  Am(A); 

where  A_(A)  is  maximized  for  A  =  2.6,  and  A  (2.6)  =  0.4871. 
in  ui 

Systems  (10)  and  (12)  —  Uniqueness 

We  will  show  that  both  systems  (10)  and  (12)  have  unique  solutions  in  the 

class 

f  |x(d)|  | 


Ej  *  <  X  :  sup  - <  ®>  (44) 

(  dtF  d2  ) 


As  in  the  case  of  systems  (7)  and  (9)  in  example  1,  if  we  show  uniqueness  for 
system  (12),  then  the  uniqueness  for  system  (10)  follows  from  theorem  A. 4, 
part  (il). 

L2 

According  to  theorem  A. 2,  the  fixed  point  S  is  unique  in  the  class 

L2  (  |x(d)| 

E  *  -  \  X  :  sup  — 7- —  <  " 

t  deF  g 

1*2  l2 
provided  that  Y°  e  E„  .  Since,  by  construction,  Y°  e  E2,  S  will  be 
1*2  r 

unique  in  E2,  If  E*  =  E2.  To  show  that  the  latter  equation  holds,  we  proceed 
as  follows. 

Let  Sn  =  L2(0),  where  L2(0)  is  the  n-th  power  of  the  operator  L2,  acting 
on  the  zero  point.  Clearly, 

L2 

s^(d)  =  b^  *  e  >  0,  for  every  deF  (45) 

Also,  it  can  be  easily  shown  by  induction  that. 


Sn(  d)=n(  (  d-A)E{N  |  A}+E{oj  |  A}+E{  ^1  A)  )-  ------ (  E{  6  |  A)  -E  { 1 1  A)  )E(n| A} 


for  every  d  e  F,  n  >1,  such  that  d  >  nA  .  EOr  d  >  2A,  letting^3' 

n  =1  -  I  -  1  in  (46),  and  using  the  fact  that  I  -  1  >  -  -  1,  yields, 

L.  A  J  *-  a  J  A 


s  (d)  > 


where  a  >  0.  (The  expressions  for  the  coefficients  u,  d,  Y  are  not  of 

interest  and,  therefore,  are  omitted) . 
l2 

If  S*  is  the  principal  solution  of  L2,  then  from  theorem  A. 1  we  have 


l2 

s*  (d)  >  sn(d)  >  0  ,  for  every  d  e  F  ,  n  *  1  (48) 

From  (45)  and  (48)  we  have  that. 


s*  (d)  e  >  0  ,  for  every  d  e  F  (49) 


From  (47)  and  (49)  we  conclude  that. 


l2 

s*  (d) 

inf  - >  0  (50) 

de  F  d2 


L2  h2 

From  (43),  and  the  fact  that  S*  <  S  ,  we  have. 


l2 

s*  (d) 

sup  - <  00  (51) 

dtF  d2 


l2 

Finally,  from  (50),  (51),  and  lemma  A.1  we  have  that 


3.  I  al  denotes  the  maximum  integer  not  exceeding  a 


i] 


3.  Stack  Algortihms 


System  ( 14 )  —  Existence  and  Initial  Bounds 


System  (14)  corresponds  to  an  operator  Li  with  bp  =  bp,  Cpv  -  cpv, 

4  >  2,  v  >2,  where  the  b^'s  and  c^v's  are  as  defined  in  (14).  Let 
X°  =  {x°(k)>  with  x°(k)  =  a^k  +  &u,  k  >  2,  and  let  Y°  =  {y°(k)>  with 
y°(k)  *  o^k  +  k  >  2.  Then,  by  straightforward  manipulations  (see  [12] 
for  details)  we  can  show  that  (20),  and  (21)  are  satisfied  (with  Lj  =  L)  if 


X  <  XQ(n)  =  ((8n3-7n2+2n+1)  -  3n+1 )/(2n(n-l) ) 


Oy  =  n( l+n+A(n-1 )  )/(2(n-1-nX)-X(n-1 )  ( 1+Xn) )  ,du  *  ( l+nXo^/tn-l ) 


otjj  *»  n/(n-1-nX)  ,  3 jj,  =  ( l+no^X)/(n-1 ) 


Thus,  if  X  <  X  (n),  then  system  (14)  has  a  solution  S  ,  such  that 


0  <  Qijjk  +  <  s  (k)  <  o^k  +  Su  ,  k  >2 


Where  o^,,  8U  are  as  given  by  (53),  and  c»£,  are  as  given  by  (54) 
We  note  that  XQ(n)  attains  a  maximum  at  n  =  3,  with  Xq(3)  =  0.3874. 


System  (14)  —  Uniqueness: 

L1  L1 

We  will  show  that  S  is  unique  in  the  class  E2^E  ,  which  is  defined 
as  follows. 


j  l*(i>|  ) 

<  X  :  sup  - <  00  ( 

(  i  >2  i2  ) 


In  order  to  do  so,  we  use  theorem  A.  4  with  Li  ^  0  1  =  ^2  anc*  =  9i(k),  where 
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g^(k)  =  k2  +  u^k  +  uQ  ,  k  >2  (56) 

and  the  coefficients  u^,  uQ  are  chosen  so  as  to  satisfy  the  conditions  stated 
in  the  theorem. 

Conditions  (a)  and  (c)  are  obviously  satisfied  for  any  uj.,  u2.  To  satisfy 
conditions  (b),  (d),  and  (e),  it  suffices  to  choose  Uj,  uQ  so  that, 

g^(i)  >  0,  and  G^(i)  >  0  ,  for  every  i  >  2  (57) 

where, 

A  “ 

G.J (i)  =  g-j(i)  -^E2CiJcg1  (k)=h(i)+ci1  +  (ci1-ni)u1  +  (ci0+ci1-n+1  )u2,i  >2  (58) 

and 

1  -  1 

h(i)  =  (1  -  -)i2  -  (1  -  2X)i  -  nA( 1+A)  ,  i  >  2 

n  n 

Given  A  <  AQ(n),  then  it  is  straightforward  to  show  that  there  exist  u^, 
u2  that  satisfy  (57).  The  following  choice  is  adopted  from  [12]: 

ui  =  max{-l,  (h2/(n-1  )+v/(nA/(n-1  )-v)  i+  1  ,  u0  =  (u-j+1)v  (59) 

where  v=  ( A(n-1 )+2)/( ( 1+A) (n-1 )+2) . 


1 


£ 


In  section  4.2,  we  will  show  that  the  sequence  of  the  conditional  mean 

L1 

session  lengths,  coincides  with  the  soluton  S  of  system 

(14),  if  this  solution  exists,  (i.e.,  if  A  <  AQ(n)).  With  this  in  mind, 

l2  l1  l2 

system  (17)  corresponds  to  an  operator  L2  with  =  b^'tS  ),  cMV  =  Cyv, 

V  A  2,  v  A  2,  where  the  b^'f-J's  and  c^'s  are  as  defined  in  (17). 

Let  X°  =  {x°(k)}  with  x°(k)  =  Yu  g-j(k),  k  A  2,  where  g^(k)  is  as 
defined  in  (56)  and  Yu  is  a  positive  real. 

Given  A  <  AQ(n),  we  choose  Yu  so  that  X°  A  L^X0)  A  0,  or  equivalently, 

Li 

Yu  g-|(i)  >b^(S  )  +  Z  cik  Yu  3l(k)  >0  ,  i  >2  (60) 

From  (55)  we  have  that  b^(S  )  <  b| l t auk+tiu^ k  >2  *  •  *■  *  2'  thus'  (60) 
holds  if,  for  every  i  A  2,  we  have, 

g-| ( i)  >  0  ,  Yu  G'j ( i)  A  b^  ({«uk+|iu>k^2)  r  i  >2  (61) 

where  G^(i)  is  as  given  by  (58). 

With  u1 ,  u2  chosen  as  in  (59),  it  follows  from  (58)  that  inequalities  (61)  are 
satisfied  if. 


bi<K*+lV>kA2> 


Yu  =  sup  -  (62) 

i  »2  G,(i) 


(it  can  be  readily  shown  that  0  <  Yu  <  “) • 

Thus,  from  lemma  A. 2,  we  have  that  if  A  <  AQ(n),  then  system  (17)  has  a 

l2 

solution  S  such  that. 


0  <  «  <i>  <  1  *  6u  1  +  4u  '  1  '  2  (63) 

where  <$u  =  u.)  Yu,  4U  =  uQ  Yu,  and  where  uj,  uQ  are  as  given  by  (59)  and  Yu 
is  given  by  (62). 

Using  similar  arguments,  one  can  readily  establish  the  following  lower 
L2 

bounds  on  S  : 

2  L2 

0  <  Yfc  i  +  <5^  i  +  Qi  <  s  (i)  ,  i  i  2  (64) 


where. 


b|({a4k+b£}k^) 


Y*  =  inf  -  ,  t>t  =  6U  ,  =  Cu  (65) 

i  32  G,  (  i  ) 


and  where  are  as  given  by  (54). 


System  (17)  —  Uniqueness 

In  thoerem  A. 3,  let  0^  =  Lj,  C>2  =L2,  and  g(i)  =  Yui  2  ♦  6U  i  +  Cu,  i  *  2, 
where  Yu,  <$u,  Cu  are  as  given  in  (63).  Then,  all  conditions  in  the  theorem 

l2 

are  satisfied.  Thus,  S  is  unique  in  the  class. 


L1  l  | x( i )  I 

Eg  =  <  X  :  sup  -  <  00 

(  Li 

s  (i)+g(i) 


L1 

It  is  clear  from  (55)  and  lemma  A.1  that  Eg  =  ^2'  where  the  class  E 2  is  as 
defined  in  the  proof  of  uniqueness  for  system  (14). 


4.2  Step  2 

In  step  1,  we  have  established  conditions  for  the  existence  of  nonnegative 


solutions  to  the  systems  of  interest,  and  we  have  identified  classes  of 
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sequences  in  which  these  solutions  are  unique.  Here,  we  show  that  the  algorithmic 
sequences  {H^},  {w^j,  where  Hj^  =  E  {h^}  and  Wi  =  E  {w^i ,  belong  to  the 
corresponding  identified  class,  and  therefore,  coincide  with  the  unique  solution 
in  the  class.  The  proof  is  based  on  theorem  A. 6,  and  is  the  same  for  all  three 
algorithms. 

Et>r  the  case  of  the  sequence  {h^},  let,  in  theorem  A. 6,  L  =  I*j, 

=  h^,  and  X1}  =  min(h^,  n),  n=  1,2,3,...  .  By  definition,  the  X^'s  and 

X"*s  satisfy  condition  (a)  in  the  theorem.  Condition  (b)  follows  from  the 

fact  that  xj?  4  n  a.e.  .  Finally,  condition  (c)  follows  from  the  operation 

of  the  algorithm.  Thus,  {H^  =  S 

l2 

Similarly,  to  show  that  IW^}  =  S  ,  we  apply  theorem  A. 6,  with  L  =  L 2, 

X^  =  w^,  and  X^  =  minCw^  n),  n  =  3,2,3,...  . 

4 . 3  Step  3 

o  o 

In  step  1 ,  we  have  already  found  upper  and  lower  bounds ,  X  and  Y  , 
respectively,  on  the  solutions  to  the  systems  of  interest.  These  bounds 
can  be  improved  either  by  computing  the  power  sequences  of  the  corresponding 
operators  with  initial  points  the  bounds  X°  and  Y°,  (lemma  A. 2),  or  by  solving 
finite  systems  of  linear  equations  that  are  truncations  of  the  original  in¬ 
finite  systems,  (theorem  A.5).  Both  methods  can  provide  arbitrarily  tight 
upper  and  lower  bounds.  We  use  the  first  method  in  the  "0.487”  algorithm, 
and  the  second  method  in  the  controlled  ALOHA  and  the  n-ary  stack  algorithm. 

1 .  Controlled  ALOHA 

FOr  system  (5),  we  apply  theorem  A.5  with  L  =  L ^  and. 


u  (i)  =  aui  +  E*u 


,  i  e  Nc 


£  (i)  =  a^i  +  3$,  ,  i  e  N0 

Rj  =  (0, 1,2. . . , j}  ,  j  E  Nq 

;where  au,  t>u,  and  a £,  8^  are  as  given  by  (27)  and  (28),  respectively.  Note 
that,  for  given  ,  Pj  is  a  finite  set  and,  therefore,  all  conditions  in  the  theorem 
are  satisfied.  Thus,  for  X  <  e-1, 

4>j  L1  Oj 

s  (i)  <  H±  =  s  ( i)  <  s  (i)  ,  0  <  i  <  j 

0.  $  . 

j  j 

where  {s  (i)}o<i<g  and  (i))o<i<j  are  the  unique  solutions  of 
the  ( j+l ) -dimensional  systems  (66)  and  (67),  respectively. 

0j  j  L1 

+  £  cik  h“  ,  i  0  <  i  <  j  (66) 

k=0 

£  j  Ll  £ 

Hi  =  bi  +  1  cik  Hi  '  i  0  <  i  <  j  (67) 

k=0 

L 2 

where  bj^  ,  bj^  are  as  defined  in  the  theorem  with  Pi  =  °i  =  b^  ,  0  <  i  <  j . 

therein.  We  solved  systems  (66)  and  (67)  for  j  =  50.  The  resulted  upper  bound 
u  £ 

Hq  and  lower  bound  HQ  on  the  mean  session  length  HQ,  can  be  found  in  table  1, 
for  different  values  of  A,  (A  <  e  ^).  For  system  (7)  we  followed  the  procedure 
described  above  with. 


L  =  L' 


u  (i)  =  (u  i  +  «u  1  +  ^u  '  1  K  No 


£  (i)  =  Yjj  i  +  Ojj  i  +  <,£  ,  i  e  Nc 


Pi  =  «i  =  bi  ,  i  e-  Pj  =  10, 1 , 2, . . . ,  j }  ,  3  e  nc 


’•  .*•  .**  »*•  •*  .*•  .  *  .  *  .** 
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where  Yu,  <SU,  £  are  as  given  by  (30),  and  Y^,  £,£  are  as  given  by  (31).  The 

resulting  bounds  W^,  on  the  mean  cumulative  dealy  Wq  are  included  in  table  1; 
they  were  computed  using  j  =  50.  Prom  the  regeneration  theorem  and  (3)  we  have^4^ 


W 


o 


D  =  -  +0.5  (68) 

AHo 

The  upper  bound  Du  =  W^/(Ah^)  +  0.5,  and  the  lower  bound  =  W^/(Ah^)+0.5 
on  D  are  included  in  table  1.  Note  that,  according  to  theorem  A. 5,  arbitrarily 
tight  bounds  can  be  obtained  by  increasing  j.  From  a  theoretical  view  point  the 
bounds  become  exact  as  j  +  °°. 


2.  The  "0.487"  Algorithm 

L1  r 

From  section  4.2  we  have  that,  for  A  <  0.487,  Hd  =  s  (d),  d  e  r,  and 
L2  H  L2 

wd  =  s  (d),  d  e  r,  where  S  and  S  are  the  fixed  points  identified  in 
section  4.1.  According  to  Lemma  A. 2  we  have  that. 


L1 

L"(Y^)  <  S  <  L^(X^)  ,  n=1 , 2, . . . ,  d  e  F  (69) 

l2 

l£(Y2°)  <  S  <  L2(x2)  '  n= 1 , 2 , • • • >  d  e  F  (70) 

where  X°  =  (c^d  +  3q}  d  £  F,  Y,°  =  {cyi  +  d  £  p 

X2  “  {Yu  d2  +  6u  d  +  V  d  t  F'  ^2  =  {^u  d2  +  6u  d  +  V  d  e  F 


and  where  au,  8U,  u^,  y^,  <5^,  ,  Yu,  6U,  are  as  given  by  (4U),  (41),  and 

(43).  Fbr  n  =  1,  and  d  =  1,  (69)  yields  the  following  bounds  on  the  mean  session 
length  H.j: 

H^|  <i  H1  <  Hf 


4  . 


The  additional  0.5  units  of  time  represent  the  mean  delay  of  a  packet,  until 
the  beginning  of  the  first  slot  following  its  arrival.  (See  footnote  1). 


where 


H,=EU|l}  +  otu(l-(l+A)e~A  +E{«.|  1>-E{5j  1}  )+iiu(l-(l+A)e  ~k) 

H*  =  H*  -  (rfu  -  3^)  ( 1- ( 1  +  A)e"A) 

The  above  bounds  can  be  found  in  table  2,  for  different  values  of  A,  (  A  <  0.487). 
for  n  =  1/  and  d  =  1,  (70)  yields  the  following  bounds  on  the  mean  cumulative  delay 
over  a  session  Wj: 

<  W1  < 

where 

W^Efal  1}+E{t|  l}+Yu(1-(1+A)e“X+E{<&-«)2|  1  > -2Et 6-A |  1} 

+  6u(l-(l+A)e-x  -E{6-A|  1>)+<;u(l-(l+A)e“X) 

wj  =  -(Cu-C£)(l-(l+A)e“A) 

u  i 

The  bounds  and  are  included  in  table  2.  From  the  regeneration  theorem 
we  have  D  =  W^/(  Ahi  ) .  The  upper  bound  Du  =  W^/(  Ah^)  ,  and  the  lower 
bound  /( Ah^)  on  the  mean  packet  delay  D  are  included  in  table  2, 

and  are  plotted  in  figure  2. 

Finally,  we  note  that  tighter  bounds  can  be  obtained  either  by  evaluating  the 
bounds  given  by  (69)  and  (70)  for  higher  values  of  n,  or  by  the  method  of  truncated 
systems  used  in  the  previous  example.  In  both  methods,  however,  we  must  first 
compute  the  conditional  probabilities  p(6,£|x)  defined  in  (9),  which  is  a 
computationally  complex  task.  Note  that  for  the  found  bounds,  (i.e.,  for  n  =  1  in 
(69)  and  (70)),  such  a  computation  is  not  required. 

3.  n-ary  Stack  Algorithm 

For  system  (14),  we  apply  theorem  A. 5  with  L  =  Li  and 


u  (i)  =  0^1+  tsu 


i  >  2 


L1 

£  (i)  =  i+  ,  i  >  2 


Pj  =  {2,3,..., j>  ,  3  >  2 


«.i  =  °i  =  b-L  ,  i  *  2 


where  c^,  and  a£,  are  as  given  by  (53)  and  (54),  respectively.  As  in  the 

®  j 

case  of  the  ALOHA  algorithm,  the  fixed  points  S  and  S  ,  in  theorem  A. 5, 

are  obtained  by  solving  two  ( j-1 ) -dimensional  systems  of  linear  equations.  If  we 
*3  0j 

denote  s  (i)  by  H^,  i  >2,  and  S  (i)  by  h“,  i  >2,  then  using  (15)  and 
(55)  we  obtain, 

H*  <  H  <  Hu 

where 

j  00 

HU  =  P0  +  P,  +  l  Pi  h£  +  E  p.  (aui  +  tSu) 
i=2  i= j+1 


H*  =  pQ  +  P,  +  l  Pi  hJ  +  E  Pi(a£  i+  S£) 
i=2  i=j+1 


VI  Jl 

In  table  3,  we  give  the  bounds  H  ,  H  for  different  values  of  A  in  the 
interval  (0,  AQ(n)),  and  for  n  =  2,3...,  where  A0(n)  is  as  given  by  (52); 
these  bounds  were  computed  using  j  =  15.  Ftor  system  (17),  we  use  theorem  A. 5 
with  L  =  L2»  and 
^2 

U  (i)  =  rui2  +  oui  +  <*u  ,  i  A  2 

1>2 

l  (i)  =  Y*i2  ♦  V  +  ,  i  >  2 

i»2 

-l>2 

l2 

=  bi<^>k*)> 


i  >  2 


where 


Yu,  6U'  ^u'  are  as  9iven  in  (62),  (63),  (65), 

l2 

(•)  =  b|(>),  with  b|(-)  as  defined  in  (17) 

Qo  =  Q?  =  If  q£  =  hJ,  2  <  k  <  j;  Q*  =  a£k  +  h£,  k  >  j+1 

qS  -  Qi  -  1;  Qfc  -  hJJ,  2  <  k  <  jf  q£  =  auk  +  eu,  k  >  j+1 

$j  Gj 

As  in  the  case  of  system  (14),  the  fixed  points  S  and  S  ,  in  theorem  A. 5,  are 
obtained  by  solving  the  corresponding  ( j— 1 ) -dimensional  systems  of  linear  equations 

Qj 

If  we  denote  s  (i)  by  W^,  i  >  2,  and  s  (i)  by  W^u,  i  >  2,  then  using  (18), 

(63),  and  (64)  we  obtain  the  following  bounds  on  the  mean  cumulative  delay  over  a 

session,  S: 

SA  <  S  <  Su 

where 

s“  -  P,  +  r  Pi  wf  +  £  Pl(Yui2  *  «ui  *  tu) 

i=2  i= j+1 

j  °° 

Sl  =  Pl  +  Z  Pi  W*  +  I  +  S^i  +  C*) 

i=2  i= j+1 

ii  £ 

The  bounds  S  ,  S  that  correspond  to  j  =  15  are  included  in  table  3.  Brora 
the  regeneration  theorem  we  have^4^  D  =  S/(Ah)  +  0.5.  The  upper  bound 
Du  =  Su/(Ah^)  +  0.5,  and  the  lower  bound  =  sV(Ahu)  +  0.5  on  the 
mean  packet  delay  D  are  included  in  table  3.  The  bounds  for  n  =  3  are  also 
plotted  in  figure  2.  From  the  table  we  see  that  bounds  found  with  j  =  15 
practically  coincide  even  for  A  close  to  AQ(n).  (According  to  theorem  A. 5, 


EE 


Finally,  note  that  the  algorithm  with  n  =  3  has  uniformly  better  mean  delay 
characteristics,  as  compared  to  the  algorithm  with  n  =  2.  From  the  operation  of  the 
algorithm  and  the  fact  that  the  quantity  XQ(n)  decreases  monotonically  for  n  >4, 
we  have  every  reason  to  believe  that  n  =  3  is  the  best  choice  among  all  n. 

5.  CONCLUSIONS  AND  PRIOR  WORK 

In  this  paper  we  have  introduced  a  method  for  the  delay  analysis  of  RMAAS,  in 
which  the  induced  packet  delay  process  is  regenerative,  and  we  have  demonstrated 
its  wide  applicability  by  applying  it  to  three  specific  examples.  The  method  is 
based  on  a  well  known  result  from  the  theory  of  regenerative  processes,  which  relates 
the  asymptotic  statistics  of  such  processes  to  quantities  that  refer  only  to  one 
cycle  of  the  process.  The  per  cycle  quantities,  (e.g.,  mean  cycle  length,  expectation 
of  the  sum  of  the  values  of  the  process  over  a  cycle) ,  are  evaluated  from  the 
solution  of  infinite  dimensional  systems  of  linear  equations.  In  Appendix  A,  we 
have  given  a  number  of  general  results  concerning  the  existence,  uniqueness  and 
approximation  of  the  solution  of  such  systems,  which  are  of  independent  interest. 

Most  of  these  results  are  generalizations  and  extensions  of  results  that  can  be 
found  in  the  early  reference  [21]. 

In  applying  the  method  to  the  three  example-algorithms,  we  have  put  emphasis 
on  the  methodology  and  rigorous  derivations  rather  than  finding  short  cuts  in  the 
analysis  of  a  particular  algorithm.  In  doing  so,  the  essential  simplicity  of  the 
method  may  have  been  obscured.  However,  to  appreciate  the  simplicity  of  the  method, 
we  note  that  only  by  using  Lemma  A. 2,  one  can  obtain  with  minimal  effort: 

1)  A  lower  bound  on  the  maximum  input  rate  that  an  algorithm  maintains  with 

finite  delay,  (i.e.,  a  lower  bound  on  the  maximum  stable  throughput  induced 
by  the  algorithm).  Note  that  for  the  first  two  examples  of  this  paper,  the 
found  bound  coincides  with  the  maximum  stable  throughput?  for  the  third 
example,  the  found  bound  is  very  close  to  the  maximum  stable  throughput 
[12],  and  since  the  induced  delays  are  already  at  very  high  values,  deter¬ 
mining  the  exact  maximum  stable  throughput  is  of  theoretical  interest  only. 


2)  Optimal  algorithmic  parameter  choices  (e.g.,  the  retransmission  proba¬ 
bility  policy  in  the  ALOHA  algorithm,  the  window  size  b  in  the  "0.487" 
algorithm,  the  splitting  parameter  n  in  the  stack  algorithm) . 

3)  Initial  bounds  on  the  mean  packet  delay,  that  can  be  used  (if  so  desired) 
to  form  finite  linear  systems,  whose  solution  can  yield  arbitrarily  tight 
bounds  on  the  mean  packet  delay. 

In  this  paper,  we  have  given  explicit  results,  only  for  the  mean  packet  delay, 
when  the  cumulative  packet  arrival  process  is  Poisson.  These  results  can  be 
readily  extended  to  the  case  of  an  arbitrary  memoryless  arrival  process,  as  long 
as  the  regenerative  character  of  the  delay  process  is  preserved.  Moreover,  the 
method  can  be  used  to  compute  higher  moments  of  the  delay  process.  This  is  due 
to  the  fact  that  if  f(x)  =  xn,  n  >  1,  in  theorem  1,  then  the  resulting  per  cycle 
quantities  are  again  related  to  the  solution  of  the  infinite  linear  systems  [12]. 
Thus,  the  results  in  Appendix  A  are  applicable. 

The  algorithms  that  served  as  examples  in  this  paper,  have  been  analyzed  in  a 
number  of  studies.  From  the  literature  on  ALOHA-type  algorithms,  we  mention  the 
work  in  [6] ,  where  the  stability  properties  of  the  version  of  the  Ctontrolled  ALOHA 
algorithm  considered  here  have  been  studied,  using  a  Markovian  model.  The  optimal 
retransmission  policy  was  derived  in  [6]  using  Pake's  lemma,  but  the  delay  analysis 
problem  was  not  addressed. 

The  delay  characteristics  of  the  "0.487"  algorithm  have  been  studied  in  [10], 
using  a  different  approach.  In  contrast  to  the  method  in  [10],  the  method  proposed 
here  does  not  require  the  computation  of  steady-state  probabilities  of  the  under¬ 
lying  Markov  chain  and,  therefore,  it  is  computationally  simpler.  Furthermore, 
since  our  approach  is  based  on  the  asymptotic  properties  of  regenerative  processes, 
it  yields  stronger  convergence  results. 

The  delay  analysis  method  of  this  paper  was  first  applied  to  the  n-ary  stack 
algorithm  in  [12].  The  analysis  in  [12]  was  stimulated  by  the  approach  taken  in 


[11],  where  the  delay  characteristics  of  the  n  =  2  algorithm  were  evaluated, 
using  the  regenerative  formulation  used  in  this  paper,  in  conjuction  with  the 
solution  of  a  functional  equation.  The  n  =  2  algorithm  has  also  been  analyzed 
in  [13],  where  a  Markovian  model  is  adopted,  in  conjuction  with  the  solution  to 
infinite  linear  systems.  We  believe  that  the  method  proposed  here  is  simpler 
than  both  of  the  above  methods.  Finally,  we  note  that  the  study  of  infinite 
linear  systems,  in  conjuction  with  the  throughput  analysis  of  the  stack  algorithm 


was  initiated  in  [9] 
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H 

o 

Ho 

w 

O 

Wu 

o 

D 

Du 

0.05 

1.00547 

1.00547 

0.05789 

0.05789 

1.651631 

1.651631 

0. 10 

1.02463 

1.02463 

0. 141334 

0.141334 

1.87936 

1.87936 

0. 15 

1.06445 

1.06445 

0.27824 

0.27824 

2.24265 

2.24265 

0.20 

1.14015 

1.14015 

0.54174 

0.54174 

2.87576 

2.87576 

0.25 

1.29097 

1.29097 

1. 17832 

1. 17832 

4.15097 

4.15097 

0.30 

1.66015 

1.66015 

3.52360 

3.52360 

7.57485 

7.57485 

0.35 

3.65080 

3.69063 

38.27660 

41.89753 

80.13219 

33.28928 

Table  1 

Delays  for  the  Cbntrolled  ALOHA 


H, 

W. 

D 

Du 

1 

1 

1 

1 

.01 

1.00025 

1.0003 

.015258 

.015258 

1.5253 

1.5255 

.05 

1.00395 

1.00474 

.08234 

.082346 

1.6348 

1.6388 

.1 

1.025 

1.030 

.  18503 

.1859 

1.796 

1.8130 

.15 

1.060 

1.061 

.3212 

.3245 

2.000 

2.040 

.2 

1.1167 

1.11367 

.5162 

.5254 

2.270 

2.352 

.25 

1.2069 

1.240 

.8243 

.8468 

2.66 

2.80 

.3 

1.356 

1.408 

1.381 

1.434 

3.270 

3.525 

.35 

1.627 

1.710 

2.6088 

2.7423 

4.358 

4.8151 

.40 

2.2279 

2.374 

6.64  38 

6.8603 

6.779 

7.670 

.45 

4.487 

4.8536 

35.012 

37.871 

16.030 

18.754 

,47 

9.110 

9.916 

163.698 

178. 178 

35.125 

41.613 

.48 

21. 175 

23. 122 

944 . 35 

1031.12 

85 . 086 

101.452 
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n  =  2 

n 

=  3 

Hu  H 

Du  D 

Hu  H 

Du  D 

o 

* 

o 

ui 

1.00565 

1.68404 

1.00619 

1.67669 

0.10 

1.02622 

1.96913 

1.02788 

1.92764 

0.15 

1.07113 

2.44552 

1.07265 

2.30451 

0.20 

1.16167 

3.33212 

1.15587 

2.91389 

0.25 

1.35801 

5.29176 

1.31346 

4.01133 

0.30 

1.92055 

11.38280 

1.65242 

6.37571 

0.35  | 

8.22892 

87.28995 

2.69508 

13.87780 

Table  3 

Delays  for  the  SA£  and  the  SA3 


cell  It  j 


cell  It  3 


cell  It  1 


the  stack 


n  cells 


a  packet 
.  arrival 


Collision 


i  +  1 


X]=0  or  1 


Idle  or  Success 


i  +  1 


Figure  1  Illustration  of  the  rules  of  the  algorithm  using  the  stack.  Xj  denotes 
the  number  of  packet  in  cell  It  j ,  at  the  beginning  of  the  i-th  slot;  N 
denotes  the  number  of  new  arrivals  during  the  i-th  slot,  and 
K,  +  K.,  +  ...  +  K  .  -  if  >  2. 


APPENDIX  A 


We  present,  in  a  generalized  format,  some  basic  results  regarding  the 

approximate  computation  of  solutions  of  infinite  dimensionality  linear  systems 

[21 J.  Let  A  be  a  denumerable  set  of  indices,  and  let  E  be  the  space  of  sequences 

X  =  tx(k)}:  A  *  R.  Given  a  set  {  c^tR,  bb  £R,i,keA  i,  let  EL  be  the  subspace 

of  E  defined  as  follows:  EL  =  (X:  £  |c^.x(k)|<  00  VitA  }.  We  define  an  operator 

keA  lk 

L  :  EL  -*•  E  as  follows:  y(i)  =  Li(X)  =  b£  +  £  c£kx(k)  1  eA  ,xeEl.  A  point 

keA 


gL  efL^  such  that. 


SL  =  L(SL) 


( A.  1  ) 


is  called  a  fixed  point  of  the  operator  L.  (A.1)  represents  an  infinite  system  of 

linear  equations  and  a  fixed  point  is  a  solution  to  this  system.  Given  an  operator 

i  a  +1  a 

L,  we  define  its  n-th  power  L  as  follows:  L  (XQ)  =  L(XQ),  L  (XQ)  =  L(L  (XQ)), 

n=1,2,...,  provided  that  XQ  eE^,  and  Ln(XQ)  eEl,  for  every  n  >  1.  The  sequence 

A  n 

{ Xn }  =  {l  ( XQ ) } ,  n=1,2,...  is  called  the  power  sequence  of  L,  with  initial 
point  XQ.  A  fixed  point  of  L  that  is  a  pointwise  limit  of  the  power  sequence  of  L, 
with  initial  point  XQ  -  0,  is  called  a  principal  fixed  point  of  L,  and  is  denoted  P> 

S*.  An  operator  0:  E^  +E  is  called  a  majorant  of  L,  iff. 


I  L  !  J 

i  c ik  '  cik 

lbil  '*  Di 


i,  k  -  A 

i  <  A 


In  this  case,  L  is  called  a  minorant  of  0.  The  notation  X  <  X’ ,  X  *  X’ ,  X,X’tE 
means  that  x(k)  <  x’(k),  (x(k)  <  x’(k)),  k  eA  .  a  point  X  eE  is  called  positive 
t nonnegative )  iff,  0  <  X  (0  *  X).  By  |x|  we  denote  the  sequence  defined  by 


|x|(k)  -  ,x(k)|,  k  t  A.  Theorems  Al,  A 2  below  are  essentially  theorems  I,  11  ^  2 
of  L  2 1 J .  They  relate  the  existence  and  uniqueness  of  a  fixed  point  of  L,  to  the 


existence  of  a  fixed  point  of  a  majorant  G  of  L. 

Theorem  A. 1  If  G  is  a  majorant  of  L,  and  0  has  a  nonnegative  fixed  point  SU, 

then  both  G  and  L  have  principal  fixed  points  S*,  S*.  Moreover,  0  <  |s*|  <  S*  4  Su 

Theorem  A. 2  If  0  is  a  majorant  of  L,  and  G  has  a  nonnegative  fixed  point  SU, 

then  the  principal  fixed  point  S*  of  L  is  unique  in  the  class  E*  CEL(  defined  as 

follows. 


E*  =  {x  eE 


sup 

ieA 


|xU>| 

s*(i) 


(1) 


Elirthermore ,  S*  is  the  pointwise  limit  of  any  power  sequence  of  L,  with  initial  point 
any  point  m  E*. 

Theorem  A. 3  below  relates  the  existence  and  uniqueness  of  a  fixed  point  of  L, 
to  the  existence  of  a  fixed  point  of  a  majorant  G  of  L,  and  it  is  a  consequence  of  thi 

theory  of  regular  systems  [21].  Its  difference  from  theorems  Al,  A2,  lies  in  the  fact 

(j  G 

that,  under  the  stated  assumptions  in  it,  we  have,  S  =  S* • 

Theorem  A.  3  If  )  is  a  majorant  of  L,  and  ..  has  a  positive  fixed  point  S°, 
such  that, 

'o 

bi 

lilt  -  ' 

1  <  A  s  (  1  ) 

then  =  S^.  Therefore,  theorem  A. 2  holds  with  S*  replaced  by  S^. 


(1)  We  adopt  the  convention: 


0  u 

-  =  1  ,  -  -  «*,  a  >  0 

0  0 


A.  1 


The  followiny  theorem  relates  the  existence  and  uniqueness  of  a  fixed  point 
of  some  operator  0 2,  to  the  existence  and  uniqueness  of  such  a  point  for  another 
operator  0 1 ,  where  the  latter  is  not  necessarily  a  majorant  of  the  former. 
Theorem  A. 4  Let  0 0 2,  be  two  operators  such  that, 

00  0 

12  U1 

(a)  cik  >  |cik|  V  i,k  tA  ,  bA  e  [0,«>,Vi  cA 


(i)  If  0-|  has  a  fixed  point  S  ,  and  there  exists  a  sequence  g;A+  r,  such  that, 

0i 

(b)  g  +  S  >0 


(c)  E  |cik  g(k)|  <  °°,Vi  e  A 
keA 

O2  0-|  0j 

(d)  | b^  |  <  (b^  +  g(i)  -  E  cik  g(k))  M,  Vi  eA  ,  for  some  M  >  0. 

keA 

Then,02  bas  a  fixed  point. 

(ii)  If  (a),  (b),  (d)  hold,  for  g  =  0,  then  the  solution  of  (?2  is  unique  in  the 

(?1 

class  E*  ,  where  E*  is  as  defined  in  Th.  A. 2. 

(iii)  If  in  addition  to  (a),  (b),  (d),  we  have  that, 

0i  0i 

bj^  +y(  i  )-  l  cikg(  x  ) 

ktA 

(e)  g  +•  S  >0,  and  inf  -  >  0, 

i>:A  0 1 

s  ( 1 )  +g( 1 ) 

0  2  01  02 

then  the  fixed  point  S  of  02  is  unique  in  the  class  Eg  d E  ,  defined  as 

follows. 

01  J  x ( i )  | 

E  -  [X<-E  :  sup  - < 

icA  0, 

s  ( 1 )  +y(  1 ) 

02  01 

S  is  the  pointwise  limit  of  any  power  sequence  of  02,  with  initial  point  in  Eg. 


Proof 


I 


0 1  0 1  0i  0i 

Part  (i):Let  Y  =  (S  +g)  M.  Since  S  =  0 -|  ( S  )  ,  we  have  that. 


01  01 

Y  /M  -  g  =  (?i  (y  /M-g) 


or 


0 1  0-i  0L  0!  0! 

y  (i)  =  Mtbi  +g( i )  -  Z  cikg(lc))  +  Z  cik  y  (i) 

keA  keA 


(A. 2) 


From  (A. 2)  and  (b),  we  see  that  the  operator  0  with  parameters, 

- 

_  n  i 

0. 


0i  0i 

b?  =  M  (bt  +g( i )-  Z  cik  g(k)),  i  eA 


keA 


c°ik  *  cik  eA  ' 


0X 

has  a  nonnegative  fixed  point  S®  ~  Y  .  Because  of  (a)  and  (d),  0  is  a  majorant 


of  02.  Prom  theorem  A.1,  we  conclude  that  02  has  a  fixed  point. 


0  _ 


'1 


Part  (ii)  s  This  follows  from  A. 2,  by  observing  that  S*  =  M  S*  ,  and,  therefore 


p0  _  r 
t*  -  t* 

Part  (iii)  :  Under  condition  (e),  theorem  A. 3  is  applicable,  and  shows  the  uniquene^ 


'1 


0 


1 


of  the  fixed  point  in  Eg  . 


The  following  lemma  is  useful  in  identifying  the  class  within  which  the  fixed 
point  of  an  operator  is  unique,  in  the  case  where  the  solution  of  the  majorant  is  not] 
exactly  known. 

Lemma  A.l  if  S, F  : A  *  r,  and, 

(a)  s,  F  are  nonnegative 


s(i) 

(b)  sup  - <  00 

icA  (j<i> 


s(  i) 

(  c )  inf - >  0  , 

IK  A  rfU> 


A.b 


|x(i)|  |x(i)j 

then  sup - <  00 ,  iff  sup  — - - <  a0,  X  c-t  ,  i.e.  the  classes 

ieA  A  <  i>  s(i) 


| x( i)  |  r  4  |x( i)  | 

E c  =  (xeA :  sup - <  °°}  and  E r  -  (xeA:  sup - <  , 

ieA  s(i)  ieA  i<(i> 


coincide. 


Proof  For  the  "if"  part  let 
|x(i) | 


sup - -  A  <  «*  >  Jx(i)|  <  A  s(i),  i  eA 

ieA  s(i) 


(A. 3) 


Because  of  (b),  we  have. 


;{i)  <  B  ^(i),  i  eA  ,  B  < 


(A. 4) 


|x(i)| 


From  (A. 3),  (A. 4),  we  conclude  that,  |x(i)  I  <  A  B  (((i),  i  eA  or,  sup  - - -  <  A  B 

ieA  (JU) 


The  proof  of  the  "only  if"  part  is  similar. 


The  lemma  below  is  used  to  establish  the  existence  of  a  fixed  point  S  of  an 

L 


operator  L,  as  well  as  upper  and  lower  bounds  on  S  .  Its  proof,  via  induction 
is  straightforward. 


Lemma  A. 2  Let  L  be  an  operator  with  nonnegative  parameters  i.e.  : 

k  „  „  .  a  ^  _  .  .  .  _  ..o  „o  ,.rL 


L  >0  i,k  t A  , 


cik 


b^  i  0,  ie^.  if  there  exist  points  Y°,  x°  eE  , 


such  that. 


(a)  Y°  <  X° 


(b)  X°  *  L(X°)  >  n 

(c)  Y°  s  L(Y°), 


then  the  power  sequence  of  L,  with  initial  points  X°  (Y°)  decreases  (increases) 


monotonically  and  pointwise,  to  a  fixed  point  SL(SL).  Furthermore,  Y°  <  SL  s  SL  <  X 


.1. 


and  S  >  0. 


It  is  generally  difficult  to  establish  tight  bounds  on  S  ,  using  the  method 


exhibited  by  lemma  A. 2.  The  following  theorem  provides  an  alternative  method  for 
the  computation  of  such  bounds. 


.6 


Theorem  A. 5  Let  L  be  an  operator  with  nonnegative  parameters: 


cLik  >  0,  i,keA  , 


bV  >0  ,  i  EA 


Let  S  be  a  nonnegative  fixed  point  of  L,  for  which  it  is  known  that 
L  <S  4  u  ,L  ,  s  ,  U  e  EL.  Let  A^CA,  A®  be  the  complement 
Df  Aj,  and  let  ♦  j  Fj/  0j  be  the  operators  with  parameters, 

<tj  Qj  /Cik  ,  i,keAj 


cik  “  cik  =  cik 


0  ,  Otherwise 


+  £  c^k  sL(k)  ,  ieA 

(  0  .  Ofche 


,  Otherwise 


•  \ 


P  i  +JAcC*  *L(1°  '  pi  <  *1'  ieA. 


,  Otherwise 


°i  +  E  -cik  uL<k)  >  «i  *  b£,  ieA 
keAj  J 

0  ,  Otherwise 


■hen,  (a)  Fj  has  a  nonnegative  fixed  point  S  3 ,  such  that. 


JU)  ,  i  eA 


0  ,  Otherwise 


3 

(b)  5j  is  a  minorant  of  F ^ ,  and  its  principal  solution  s  *  is  such  that 


0 .  F  F 

J  j  i 
o  *  s*  <  s*  4  s 


Di 

(c)  Oj  is  a  majorant  of  Fj,  and  if  sup  - <  <*>,  then  Oj  has  a 

ieA.  L 


_  1£  J  L: 

Oj  bi 

nonnegative  fixed  point  S  ,  such  that, 

Fj  Oj  Oj 

o  <  s*  4  s*  4  s 


A.  7 


(d)  If  in  addition  to  the  previous  conditions,  also  SL  >  0,  and 


*■  r 

inf  - >  0,  then  the  operators  ♦  r  :  y-  have  respective 

ieA  sL(i)  3  3  3 

*j  Fj  °j 

unique  fixed  points,  S  ,  S  ,  S  ,  in  the  class 


r  J  (  r  | 

E  =  <  X  t.E  :  sup - <  00  >  , 

(  ieE  FJ  ) 


Fj  Gj 

and  S  <  S  t  S 


s  ( i) 


Remark 


If  P^  is  a  finite  set  with  >  0,  Vi  £  P^ ,  conditions  (c)  and  (d)  are 


clearly  satisfied.  If  in  addition,  =  cl  =  b^,  and  A^/A,  then  it  can  be 

Fj 

shown  [21]  that,  S*  +  SL,  and  S*  ♦  SL,  pointwise. 


The  quantities  of  interest  in  the  various  random  access  algorithms  are  statist! 


of  random  variables,  where  many  of  those  statistics  are  fixed  points  of  some  operatoj 
Theorem  A. 6  is  used  to  justify  the  latter  statement  and  appeared  in  [15]. 

Theorem  A. 6  Let  L  be  an  operator  with  nonnegative  parameters,  that  has  a  unique 

L  L  lxil 

nonnegative  fixed  point  S  in  the  class  Eg  =  (xeP  :  sup  -  <  "} . 

itP  g(i) 

Let  Ix^11},  Ix^},  i  e  P,  n  e  n,  be  families  of  random  variables,  such  that, 

(a)  0  <  a.e.  for  every  i  ep 

(b)  x£  4  M^gti),  a.e.  for  every  i  P,  <  00 

(c)  fn  <  L(  f  n )  ,  f  =  L(f),  where  fn(r)  =  E  Ix^),  f(i)  =  E  (x^ 

Then,  f  coincides  with  the  unique  fixed  point  SL  in  E^. 

Proof 

We  observe  that  because  of  (b),  then  f  tEL,  and  because  of  (c)  and  lemma  A. 2, 
then  fR  <  SL.  Since  also  f  increases  to  f  pointwise,  because  of  (a)  and  the 
monotone  convergence  theorem,  we  conclude  that  f  <  SL,  and  therefore  f  t 


The  assertion  now  follows,  from  the  fact  that  f  is  a  fixed  point  of  L. 


B. 


APPENDIX  B 


In  section  4.2,  we  saw  that  the  computation  of  conditional  expectations, 
E{x|d},  is  required.  In  this  appendix,  we  show  that  those  conditional  expectations 
can  be  computed  with  high  accuracy.  Let  us  define, 

E{x|d,  k}  :  The  conditional  expectation  of  the  random  variable  X,  given  that 
the  arrival  interval  contains  k  packets,  and  has  length  d. 

Then, 


E{X|d}  ■  Z  E{x|d,k>e_Xd 
k=0 


(B.l) 


Using  the  rules  of  the  algorithm,  the  quantities  E{X|d,k}  can  be  computed 
recursively  ,  as  follows. 


EU/d,k}  «=  EU/l,k};  V  deF 
EU/1,0}  -  EU/1,1}  «  1 


E{m,k} 
where  Pk 


.  k-x  . 

(l+PJ+EU/l,k-l}Pk+  Z  E{jt/l,i}p£)/(l-2Pk);  k  _>  2 

i-2 

( i ) 2_k 


(B.2) 


E{S/d,k}  -  d  E{6/l,k};  V  deF 

E{6/1,0}  -  E{6/l,k}  -  1  (B.3) 

E{6/l,k  -  (Pk+pk+E{6/l,k}P^+  Z  E{6/l,i>ph/(2-(l-Po));  k  >  2 

o  o  1=2 

EU2/d,k}«  EUz/l,k};V  deF 
EH2/!,!)}-  E{Jl2/l,l}  -  1 

E{42/l,k}-  (2E{Jl/l,k>+2E{Jl/k,k-l}*P^+P^-l+P^E{22/l,k-l}+  Z  E{E2/l,i}pJ)/ 

i-2 

(l-2pg);k>2  (B.4) 

E{<S2/d,k}-  d2F{62/l,k}  ;  V  deF 
E{62/1,0}«  E{ 62/1,1}  -  1 

E{62/l,k}-  (.25(P^+Pq)+.5  P^E{6/l,k}+.5  P*  E{6/l,k-l}+.25  E{52/l,k-l}P^+ 

k-i  .  . 

.25  Z  E{62/l,i}Pk)/(l-.5P3)  ;k>2 


(B.5) 


E{6i/d,k} 

-  d  E{6i/l,k)  ;  V  deF 

E{54/1,0} 

-  E{6£/1,1}  =  1 

E{«i/l,k} 

-  (E{6/l,k}+.5P^E{A/l,k-l}+.5P^E{i/l,k-l}+.5 

pJ+.5E{6/l,k-l>pJ+ 

+.5P^E{5)l/l,k-l}+.5k£1E{6Jl/l,i}pJ;)/(l-ph  ;  k 
A  i«2  1  0 

>  2  (B.6) 

E{N/d,k} 

i  Am 

-  E{N/l,k>;  V  deF 

e{n/i,o} 

-  0,  e{n/i,i}  -  1 

E{N/l,k} 

-  p]f-HP^E{  N/l  ,k-l}  +k£1E  {N/l , i  }P^ ;  k  >  2 

A  A  1  =2  1 

(B.7) 

E{a)/d,k} 

1 

-  E{w/l,k};  V  deF 

E{u)/1,0} 

-  0,  E{w/l,l}  *  1 

(B.8) 

E{(t)/l,k} 

k-1 

-  (p!f+E{  N/l , h}+P^E { N/l , k-l}+P^E {w/ 1 , k-1 }+  £  E{a)/l,k}pk)/(l-2p)S  ;  k  >  2 
l  x  x  i„2  1  u 

E{tj>/d,k} 

-  d  E{i|»/l,k>;  V  deF 

E{^/1,0} 

-  0,  E{ij»/l,l}  -  | 

(B.9) 

Ety/l,k} 

k-1 

-  (E{N/l,k}(l-Po)-P^E{N/l,k-l)+.5P^4E{l|>/l,k-l}P^+  £  Bty/l,i}ph 

X  X  L  i«2 

/ (2  (l-pp>  ;  k  >  2 

From  formulas  (B.2)-(B.9),  we  see  that  a  finite  number. 

M,  of  terms  from  the 

Infinite  series  (B.l) ,  can  be  easily  computed.  Also,  for  large  k  values,  and  based 
on  the  recursive  expressions,  simple  upper  and  lower  bounds  on  E{X/d,k}  can  be 

so  (Xd) 

developed.  Those  bounds  can  be  used  to  tightly  bound  the  sum  £  E{X/d,k}e  Vry 

k>*H-l 

Remark  It  can  be  also  proved  that 
E{N/d}  -  XE{6/d} 

Etip/d}  -  Xd  E{6/d}  -  XE{«2/d} 


END 
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